Nucleon Form Factors in QCD 
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We calculate the electromagnetic and the axial form factors of the nucleon within the framework 
of light cone sum rules (LCSR) to leading order in QCD and including higher twist corrections. In 
particular we motivate a certain choice for the interpolating nucleon field. We find that a simple 
model of the nucleon distribution amplitudes which deviate from their asymptotic shape, but much 
less compared to the QCD sum rule estimates, allows one to describe the data remarkably well. 
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I. INTRODUCTION 

Form factors play an extremely important role 
in the studies of the internal structure of compos- 
ite particles as the measure of charge and current 
distributions. In particular, the pioneering study 
of the nucleon form factors by Hofstadtcr and col- 
laborators 01 demonstrated that the nucleons have 
a finite size of the order of a fermi. The behavior 
of the form factors at large momentum transfers 
is especially interesting. Already in the pre-QCD 
times it vifas established that, if one can treat the 
hadrons at high momentum transfer as coUinear 
beams of N valence quarks located at small trans- 
verse separations and exchanging intermediate glu- 
ing particles with which they interact via a dimcn- 
sionless coupling constant, then the spin- averaged 
form factor behaves asymptotically as 1/{Q^)^~^ 
[3. This hard-exchange picture and the resulting 
dimensional power counting rules [^ y] can be for- 
mally extended onto other hard exclusive processes. 

After the advent of quantum chromodynamics, 
this hard-gluon-exchange picture was formalized 
with the help of the QCD factorization approach 
to exclusive processes |j, la p ■ This approach in- 
troduces the concept of hadron distribution am- 
plitudes (DAs). They arc fundamental nonpcr- 
turbativc functions describing the momentum dis- 
tributions within rare parton configurations when 
the hadron is represented by a fixed number of 
Fock constituents (quarks, antiquarks and gluons). 
It was shown that in the Q^ -^ cx) limit, form 
factors can be written in a factorizcd form, as a 
convolution of distribution amplitudes related to 
hadrons in the initial and final state times a "short- 
distance" coefficient function that is calculable in 
QCD perturbation theory. The leading contribu- 
tion corresponds to DAs with minimal possible 
number of constituents — three for baryons and 
two for mesons. 

The essential requirement for the applicability of 
this approach is a high virtuality of the exchanged 



gluons and also of the quarks inside the short dis- 
tance subprocess. More generally, in the case of the 
nucleon form factors the hard perturbative QCD 
(pQCD) contribution is only the third term of the 
factorization expansion. Schematically, one can en- 
visage the expansion of, say, the Dirac electromag- 
netic nucleon form factor Fi{Q^) of the form 



FiiQ'' 



A{Q'') 



a,(g2)\ B(Q2) 



a^my C 



(1) 



where C is a constant determined by the nucleon 
DAs, while A{Q'^) and B{Q^) arc form-factor-typc 
functions generated by contributions of low vir- 
tualities, see Fig. Q] The soft functions A{Q'^) 
and B{Q^) are purely nonperturbative and can- 
not be further simplified e.g. factorizcd in terms 
of DAs. In the light-cone formalism, they are 
determined by overlap integrals of the soft parts 
of hadronic wave functions corresponding to large 
transverse separations. Various estimates suggest 

that A(Q2) < 1/Q6^ 5(Q2) < 2/Q4 ^^^ j^^ ^j,^y 

large Q^ they are further suppressed by the Su- 
dakov form factor. To be precise, in higher orders 
in as{Q) there exist double- logarithmic contribu- 
tions ^ 1/Q^ [3 that arc not factorizcd in the stan- 
dard manner; however, also they are suppressed by 
the Sudakov mechanism [alsllld- Thus, the third 
term in ^ is formally the leading one at large Q^ 
to power accuracy. 

The main problem of the pQCD approach is a 
numerical suppression of each hard gluon exchange 
by the as/7r factor which is a standard perturba- 
tion theory penalty for each extra loop. If, say, 
as/TT ~ 0.1, the pQCD contribution to baryon 
form factors is suppressed by a factor of 100 com- 
pared to the purely soft term. As the result, the 
onset of the perturbative regime is postponed to 
very large momentum transfers since the factoriz- 
able pQCD contribution 0{l/Q^) has to win over 
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FIG. 1: Structure of QCD factorization for baryon form factors. 



nonpcrturbative effects that are suppressed by ex- 
tra powers of 1/Q^, but do not involve small coef- 
ficients. There is a growing consensus that "soft" 
contributions play the dominant role at present en- 
ergies. Indeed, it is known for a long time that 
the use of QCD-motivated models for the wave 
functions allows one to obtain, without much ef- 
fort, soft contributions comparable in size to ex- 
perimentally observed values (see, e.g. [ill Ha . Il3| ) . 
A modern trend [ij, [iSj is to use the concept 
of generalized parton distributions (GPDs) to de- 
scribe/parametrize soft contributions in various ex- 
clusive reactions, see [l6|, [iTl ll8| for recent re- 
views, and the models of GPDs usually are chosen 
such that the experimental data on form factors 
are described by the soft contributions alone, cf. 
Refs. [Ulil ili. A subtle point for these semi- 
phenomenological approaches is to avoid double 
counting of hard rescattering contributions "hid- 
den" in the model-dependent hadron wave func- 
tions or GPD parametrizations. 

The dominant role of the soft contribution for 
the pion form factor at moderate momentum trans- 
fers, up to Q^ ~ 2 — 3 GeV^, is supported by 
its calculation \22L l23l | within the QCD sum rule 
approach |2J|. The application of the method at 
higher Q^ faces the problem that the inclusion of 
nonpcrturbative effects due to vacuum condensates 
through the expansion over inverse powers of the 
Borel parameter M^ interferes with the large-Q^ 
expansion of the form factors, producing an ill- 
behaved series of the type J2n '^niQ^ / ^^^)"' ■ For 
the nucleon form factors, the QCD sum rule ap- 
proach only works in the region of small momen- 
tum transfers Q'^ < 1 GeV^ [H H^. To extend 
the results to higher Q^, it was proposed [23 to re- 
sum the (Q^/M^)" contributions originating from 
the Taylor expansion of simple models for nonlocal 
condensates. Another approach _23^, ^] is to use 
the so-called local quark-hadron duality approxi- 
mation, in which the condensates are effectively 
neglected. The parameter-free results for the pion 
and nucleon form factors obtained in this way are 
in a rather good agreement with the existing data. 



We also have to mention the dispersion approach 
of Ref . [29I [3OJ which allows to analyze form factors 
for all momenta (space- and time-like) in a largely 
model-independent manner in terms of spectral 
functions on a hadronic level. Also, in future, 
one expects that the rapid development of lattice 
QCD will allow one to calculate baryon form fac- 
tors to sufficient prec ision from first principles, see 
e.g. [3^, [33, [3^ [33. Such studies are necessary 
and interesting in its own right, but do not add to 
our understanding of how QCD actually "works" 
to transfer the large momentum along the nucleon 
constituents, the quarks and gluons. The main mo- 
tivation to study "hard" processes has always been 
to understand hadron properties in terms of quark 
and giuon degrees of freedom; for example, the ra- 
tionale for the continuing measurements of the to- 
tal inclusive cross section in deep inelastic scatter- 
ing is to extract quark and gluon parton distribu- 
tions. Similar, experimental measurements of the 
form factors at large momentum transfers should 
eventually allow one to determine baryon distribu- 
tion amplitudes and this task is obscured by the 
presence of large "soft" contributions which have 
to be subtracted. 

In Ref. [35J we have suggested to calculate 
baryon form factors for moderatel y la rge Q^ using 
light-cone sum rules (LCSR) [3ll33- This tech- 
nique is attractive because in LCSRs "soft" contri- 
butions to the form factors are calculated in terms 
of the same DAs that enter the pQCD calculation 
and there is no double counting. Thus, the LCSRs 
provide one with the most direct relation of the 
hadron form factors and distribution amplitudes 
that is available at present, with no other nonpcr- 
turbative parameters. 

The basic object of the LCSR approach is the 
correlation function 



dxe-'^«-(0|r{7y(0)j(x)}|iV(P)) 



in which j represents the electromagnetic (or weak) 
probe and rj is a suitable operator with nucleon 
quantum numbers. The other (in this example, 
initial state) nucleon is explicitly represented by its 



state vector \N{P)), see a schematic representation 
in Fig. [3 When both the momentum transfer Q^ 




FIG. 2: Schematic structure of the light-cone sum rule 
for baryon form factors. 

and the momentum (P')^ = {P — q)^ flowing in 
the T] vertex are large and negative, the asymp- 
totics of the correlation function is governed by 
the light-cone kinematics x^ — *■ and can be stud- 
ied using the operator product expansion (OPE) 
T{77(0)j(x)} ~ Y.Ci{x)0^{{)) on the light-cone 
x^ =■ 0. The x^-singularity of a particular per- 
turbatively calculable short-distance factor Ci (x) is 
determined by the twist of the relevant composite 
operator d, whose matrix element (0|C'i(0)|A^(P)) 
is given by an appropriate moment of the nucleon 
DA. Next, one can represent the answer in form 
of the dispersion integral in (P')^ and define the 
nucleon contribution by the cutoff in the quark- 
antiquark invariant mass, the so-called interval of 
duality sq (or continuum threshold). The main role 
of the interval of duality is that it does not allow 
large momenta |fc^| > sq to flow through the rj- 
vertex; to the lowest order 0{a^) one obtains a 
purely soft contribution to the form factor as a 
sum of terms ordered by twist of the relevant op- 
erators and hence including both the leading- and 
the higher- twist nucleon DAs. Note that, in dif- 
ference to the hard mechanism, the contribution of 
higher-twist DAs is only suppressed by powers of 
|(P')^I ~ 1 - 2 GeV^ (which is translated to the 
suppression by powers of the Borel parameter after 
applying the usual QCD sum rule machinery), but 
not by powers of Q^ . This feature is in agreement 
with the common wisdom that soft contributions 
are not constrained to small transverse separations. 
The LCSR expansion also contains terms gener- 
ating the asymptotic pQCD contributions. They 
appear at proper order in a^, i.e., in the 0{as) 
term for the pion form factor, at the a^ order for 
the nucleon form factors, etc. In the pion case, 
it was explicitly demonstrated [33, |39j that the 
contribution of hard rescattering is correctly repro- 
duced in the LCSR approach as a part of the 0{as) 



correction. It should be noted that the diagrams 
of LCSR that contain the "hard" pQCD contri- 
butions also possess "soft" parts, i.e., one should 
perform a separation of "hard" and "soft" terms 
inside each diagram. As a result, the distinction 
between "hard" and "soft" contributions appears 
to be scale- and scheme-dependent [33 ■ During the 
last years there have been numerous applications of 
LCSRs to mesons, see |4y,|4l| for a review. Follow- 
ing the work J35j nucleon form factors were further 
considered in this framework in Refs. [42, |43, |4^ 
and the weak decay A^ —^ piui in [43 ■ The gen- 
eralization to the A''7A transition form factor was 
worked out in |4y| . 

In this paper, we go beyond the original work 
[33 in several important aspects. First, we present 
a detailed study using different interpolating cur- 
rents for the nucleon and choose one which ap- 
pears to be the optimal. Second, we calculate both 
the electromagnetic and weak decay form factors. 
Third, we make an update of the parameters of 
higher-twist DAs which feature prominently in this 
approach and are important numerically. We then 
formulate a simple model for the DAs that pro- 
vides a good description of the available experimen- 
tal data. Finally, we include a complete summary 
of higher-twist DAs and work out the light-cone 
expansion of three-quark currents for all Lorentz 
structures, which extends the results given in [33 ■ 
The presentation is organized as follows. Section 
2 is introductory and contains the form factor def- 
initions and some general discussion. Section 3 is 
devoted to the construction of the sum rules. We 
present here our numerical results and the conclu- 
sions. The final section 4 is reserved for a short 
summary and an outlook. The paper contains five 
appendices devoted to a summary of correlation 
functions to tree-level accuracy, summary of three- 
quark nucleon distribution amplitudes, the OPE of 
generic three-quark amplitudes to twist-5 accuracy, 
and the QCD sum rule estimates for the parame- 
ters of higher-twist DAs. 



II. PRELIMINARIES 

A. Electromagnetic Form Factors of the 
Nucleon 

The matrix element of the electromagnetic cur- 
rent 

i""(-^) = e„i7(a;)7pw(x) -|- edJ(x)7^d(x) (2) 

taken between nucleon states is conventionally 
written in terms of the Dirac and Pauli form factors 



Fi(Q2) and^2(Q'): 

{N{P')\f;^\0)\NiP)) = 

^2^ ■^ti.i'ft 



-^N{P') 



7t.Fi{Q' 



2m 



N 



F2m 



7V(P),(3) 



where P^ is the initial nucleoli momentum, P^ = 



i-jV 



, P' = P - q, Q^ := -q^, ct^^ = §[7^,7.^] and 
N{P) is the nucleon spinor. 

Experimental data on the scattering of electrons 
off nucleons, e.g. e^ + p ^ e^ + p, is usually pre- 
sented in terms of the electric Ge{Q^) and mag- 
netic Gm{Q'^) Sachs form factors which are related 
to Fi^2{Q^) as 



GA/(g2)=Fi(Q2) + ^2(g2)^ 



GB(g2)=Fi(Q2) 



4m% 



F2m- 



(4) 
(5) 



In the Breit frame the form factors Ge{Q^) and 
Gm{Q^) can be thought of, loosely speaking, as 
the Fourier transforms of the charge distribution 
and magnetization density in the nucleon. 

The normalization of the form factors at Q^ = 
is given by the nucleon charges and magnetic mo- 
ments 

G^O) = 1, Glj{Q) = fip = 2.792847337(29) , 

Gb(0) = 0, Gm(0) = Hn = -1.91304272(45) (6) 

for the proton and the neut ron, res pect ively | 4 7| . 
Experimentally one finds 0, ligllsol. I5l ll52ll53l 

|5^;5^,5S|53,|53| that the magnetic form factors of 
both the proton and the neutron can be described 
by the famous dipole formula: 



GlriQ')-' 



Gm{Q 



^J■p 






{^-%y 



(7) 
(8) 



with fj,Q — 0.71 GcV^. For the electric form fac- 
tor of the neutron the measured values are close 
to zero |59|, |60j. In the case of the electric form 
factor of the proton the situation was controver- 
sial for some time, with the experimental mea- 
surements using the classical method of Rosen- 
bluth separation producing very different results 
compared to the ones obtained using the method 
of p olarization transfer. Recently it was argued 

mEZ. 3^ 3^ ^ iS(k £Ii, M, M that the for- 
mer approach is not applicable for sufficiently large 
momentum transfers, as the contribution of the 



electric form factor to the spin-averaged cross sec- 
tion is strongly contaminated by contributions of 
the two- pho ton exc hang e. The existing estimates 
IHillHlililiailEiEi indicate that 
the two-photon exchange corrections have the right 
sign and order of magnitude to bring the values ob- 
tained via the Rosenbluth separation to the ones 
measured using the polarization transfer, although 
the situation is not finally settled. In this work we 
rely on the polarization transfer data [73, IZll, lZ3 ■ 



B. Charged Weak Form Factors 

In order to describe charged current (CC) neu- 
trino reactions like 



i^fi + n^fi + p, 


(9) 


^M + P^M^ + n, 


(10) 



one has to deal with matrix elements between nu- 
cleon states of the vector V^p*-^ and the axial- vector 
current ^P*-^: 

V;C''^(x)=^7^T+*(x), (11) 

^;:'^(a:;)=%^75r+*(x), (12) 

where VP is an (iso)spinor consisting of an up- and 
a down-quark and t+ = l/2(ri + zt2) is a linear 
combination of the familiar Pauli matrices. 

One defines the corresponding vector and the 
axial- vector form factors as 

rCC 



(iv(p')inr^(o)|A^(p)) = 



=iV(P') 



7,Fr(0^)-z 



i^^i^g CC 



2mN 



F^^iQ' 



N{P), 



{N{P')\A';:^{0)\N{P))^N{P') i,G^/^{Q^) 



'^^ -G'^'^iQ') - ^5^Gg^(Q2) ^^^,(p) . 



2mN 



2771 AT 



(13) 



The vector form factors can be related with the 
electromagnetic ones with the help of isospin sym- 
metry, to wit 

{p\urd\n) = {p\jl^Jp) ~ {n\ji^Jn), 

(p|u7''75d|n) = (p|(u7''75U- d'y^"f5d)\p). (14) 

The first relation gives, e.g. 

Ff^iQ')=FfiQ')~F^{Q'), 

F^^{Q^)=F^{Q^)-Fi\Q^). (15) 



The axial form factor G'^'~'' {Q^) can be deter- 
mined either from quasi-elastic neutrino scattering 
or from pion electroproduction (with the help of 
current algebra). The neutrino data are available 
for Q2 values up to 3 GeV^ JH El El ■ They 
were reanalyzed recently in |73 ■ The pion electro- 
production data exist for Q^ < 1 GeV^, e.g. |7a |. 
After the inclusion of a finite pion mass correction 
in the analysis |79| , the extracted form factor agrees 
well with the determinations in neutrino scattering. 
All the existing experimental data for G'^'-^(Q^) at 
Q^ < 1 GeV^ are very well described by the dipole 
formula: 



Gi'^m 



9A 



i'-m 



(16) 



with gA = 1.267 ± 0.004. The mass parameter is 
fitted to be Ma = 1.001 ± 0.020 GeV and Ma = 
1.013 ± 0.015 GeV from neutrino scattering and 
pion electroproduction, respectively. 

The pseudoscalar form factor Gp^{Q^) can be 
extracted separately from muon capture of the pro- 
ton fi~ -\- p —f v^ + n OT from pion electroproduc- 
tion. In this case only the data for Q^-valucs below 
0.2 GeV^ exist |80j, which is to low for the applica- 
tion of our method. Using PCAC and the pion pole 
dominance model one can express the pseudoscalar 
form factor Gp'~'' in terms of the axial form factor 
Ua ■ 



r<cc 



iQ' 



Q2 + ^2 



(17) 



This form is consistent with the conservation of the 
flavor nonsinglct axial current in the chiral limit 

("4 ^ 0). 

Finally, the tensor form factor Gp^{Q^) must 
vanish by virtue of the isospin symmetry and T- 
invariance, so it is normally not included. The 
reason why we leave it in Eq. H13() is that in our 
approach the initial and the final state nucleons 
are treated differently, so that T-invariance is not 
manifest. 



C. Neutral Weak Form Factors 

The cross section for clastic neutrino-proton and 
neutrino-neutron scattering can be expressed in 
terms of matrix elements of a vector V,!^'~^ and an 
axial- vector A^ neutral currents: 



C(-) = 



w I u'y^u[x) 



iiin^Ow d-i^d{x) , (18) 



K''^^)=2 N'^75W(2:) - dn^dix)] , (19) 

where 9w is the Weinberg angle. The matrix el- 
ements of neutral currents over the nucleon states 
are conventionally written as 

rNC 



(iV(P')lC(0)|A^(n) 



--N{P') 



J,Fr{Q')~^^£^Fi^^{Q'^ 



2mN 



N{P) 



{NiP')\A^^iO)\N{P))^N{P') 



l.GTiQ') 



2m 



N 



criQ') 



;(^liv(l ^NC 



2mN 



G^^iQ') i,N{P) 



(20) 



The vector form factors F-l^'~^ and F^^ are, again, 
just linear combinations of the electromagnetic 
form factors of the nucleon. For the axial neu- 
tral weak form factors {Ga'^ and Gp'^) there is 
little data, and only for Q'^ < 1 GeV^ [sj which 
is below the region we try to describe theoretically. 
The tensor form factor G^'^{Q^) must vanish by 
virtue of T-invariance; the reason we include it will 
become clear later. 



D. Light-Cone Kinematics 

Having in mind the practical construction of 
light-cone sum rules that involve nucleon DAs, we 
define a light-like vector z^ by the condition 

q-z = 0, z^ = (21) 

and introduce the second light-like vector vector 



Pm=^m 



1 m% 



2 "P-z 



p2 = , (22) 



so that P ^ pii the nucleon mass can be neglected, 
niN -^ 0. The photon momentum can be written 
as 



9^ = l^t^ + ^p 



P-q 

P- z 



(23) 



We also need the projector onto the directions or- 
thogonal to p and z, 

atv = 9t^u {p^,Zu+PuZ^,), (24) 

pz 



and use the notation 



^/^^ •} 



(!•, tii.Z 



2p CiflP 7 



(25) 



for arbitrary Lorentz vectors a^. In turn, aj_ de- 
notes the generic component of a^ orthogonal to z 
and p, in particular 



qi 



% 



pq 

pz 



(26) 



We use the standard Bjorken-DrcU convention 
[821 for the metric and the Dirac matrices; in par- 
ticular, 75 = i7°7^7^7'^, and the Levi-Civita tensor 
£fj,uXa is defined as the totally antisymmetric tensor 

with £0123 = 1- 

Assume for a moment that the nucleon moves 
in the positive e^ direction, then p"*" and z~ are 
the only nonvanishing components of p and z, re- 
spectively. The infinite momentum frame can be 
visualized as the limit p^ ^ Q —>■ oo with fixed 
P ■ z = p ■ z ^ 1 where Q is the large scale in the 
process. Expanding the matrix element in powers 
of l/p"*" introduces the power counting in Q. In this 
language, twist counts the suppression in powers of 
p"*" . Similarly, the nucleon spinor Nj (P, A) has to 
be decomposed in "large" and "small" components 
as 

N^{P, A)=-^ W+ ^^) N^{P, A) 
2p ■ z 

=N+{P,\)+N~{P,X), (27) 

where we have introduced two projection operators 



A+ 



2p- z 



A- 



2p- z 



(28) 



that project onto the "plus" and "minus" compo- 
nents of the spinor. Note the useful relations 



^ N{P) = mNN+ (P) , /iN{P) 



2p- 



niN 



-N-iP) 
(29) 



that arc a consequence of the Dirac equation 
fN{P) = ■mNN{P). Using the explicit expressions 
for N{P) it is easy to sec that A+iV = A^+ -■ y/p+ 
while A^iV = iV" - l/^/p^- 

Note that all expressions are invariant under the 
reparametrization z^ — s- az^ where a is a real num- 
ber; we will use this freedom to set z^ equal to the 
"minus" component of the distance between the 
currents in the operator product. 



III. LIGHT-CONE SUM RULES FOR 
BARYON FORM FACTORS 

A. Choice of the current 

As already mentioned, in the LCSR approach 
one of the participating nuclcons is replaced by a 



suitable interpolating current for which there are 
several choices. Altogether, there exist three local 
operators with isospin / = 1/2 numbers that do 
not involve derivatives |8J| . They can be chosen as 

iu{x)=e'^''[u\x)C^^u^{x)\ i^rd\x), (30) 
772(2;) =£'■'"'= [u\x)Ca^,u^{x)\ 75a^''d'-(^),(31) 



m{x) = -e''\[u\x)C)fu^{x)\ 75/d'(a;) 



[u\x)C^d^{x)\-i^^u''{x) 



(32) 



where u{x) and d{x) are the w-quark and the d- 
quark field operators, respectively, «,j, fc are color 
indices, C is the charge conjugation matrix |82l | and 
z is a light-like vector, z^ = 0. The corresponding 
couplings 

(0| 7,1(0) |Ar(P))=AimArAr(P), 
(0| 772(0) |Ar(P))=A27nAriV(P), 

(0| 7/3(0) |iV(P))=/N(Pz) ^N{P) (33) 

are well known, albeit with limited precision, from 
the vast QCD sum rule literature, see Eq. (|B.13p 
in Appendix B for the current estimates. 

The operator 771 is known as the loffe-current 
851 . There is overwhelming evidence (see, however 
86j) that this current gives rise to more accurate 
and reliable sum rules compared to 772 [83. By 
this reason we do not use 772 in the construction 
of the LCSRs below. The both constants Ai and 
A2 do appear in the sum rules, however, as they 
determine the normalization of the highcr-twist-4 
nucleon DAs. 

In turn, the operator 773 is twist-3 and the cor- 
responding coupling /tv determines the normaliza- 
tion of the leading twist nucleon DA |83|. Using 
the currents of lower twist is in general advanta- 
geous, as the corresponding correlation functions 
have lower dimension and are less affected by the 
model-dependent continuum subtraction. LCSRs 
obtained with this current are most close in spirit 
to pQCD factorization. The price to pay is that 
this current couples both to the spin J = 1/2 and 
spin J = 3/2 baryons and it is unclear whether the 
unwanted / = 3/2 contributions to the correlation 
function are sufficiently suppressed. This current 
has been used more rarely in the practice of QCD 
sum rule calculations, so that there is less experi- 
ence. 

In Ref. |35j a modification of the current (|32|1 
was used 

774(2;) = £*■'■'= [u'{x)C^u'{x)\ 75^d'^(x) , 

(0| 7/4(0) |A^(P)) = /n(Pz) ^N{P) , (34) 



which, in difference to 773, also couples to the 
isospin / = 3/2 states (e.g. the A-isobar). 

A priory, it is not obvious whether using the 
pure isospin- 1/2 or, similarly, pure spin-f/2 cur- 
rents improves the accuracy of the sum rules. 
It is conceivable that the summation over quan- 
tum numbers makes the duality approximation for 
taking into account contributions of heavy reso- 
nances and the continuum more accurate and also 
suppresses poorly known contributions of higher- 
dimension (twist) operators. In a different context, 
there have been various proposals to add together 
correlation functions of opposite parity, |89j, use 
chirally projected quark fields, e.g. |93, etc. 

We believe that there is no general recipe; one 
has to consider each case separately and the con- 
clusions can vary. For the problem at hand, it was 
noticed in [4J| that using the currents 7^3 and 774 
one obtains numerical results for the nucleon form 
factors that differ significantly from one another. 
In this work we demonstrate that this difference 
is due to the contamination of the sum rules [35J 
by the contributions of isospin / = 3/2 states and, 
therefore, the use of the pure-isospin current r/3 is 
strongly preferred compared to 774. On the other 
hand, the sum rules obtained using 771 and 773 are 
complementary to a large extent and both of them 
arc useful. Nevertheless, using the loffe current rji 
produces the sum rules that arc more stable and 
seem to be superior in all respects, which makes 
this current to be our final choice. 



B. Correlation functions 

We consider the set of correlation functions 
r^-<=-(P, q)=t fd^xe^^^ (0| T hmn^)] \N{P)) , 



Tl''^{P,q)^t d^xe^'^^ {0\T[rn{Om{x)] \N{P)) , 



n^'{P,q)^t d^xe^^- {0\TMO)j:{x)] \N{P)) 



(35) 

where T denotes time-ordering, \N{P)) is the pro- 
ton state with four- momentum P^, P^ — mj^, rji 
with i = 1,3,4 are the nucleon currents defined in 



(|5n|) . Further, j°™ is the electromagnetic current 
defined ^ whereas j^ and j^ are the isospin-one 
axial and vector currents, respectively; 



Jt''"'ix) = -[u{x)-f^-/5u{x) - d{x)-f^-/5d{x)] , 

Jl''"'{x)=u{x)j^'-f5d{x) , 
Jl''"'{x)=-^[u{x)ji,u{x) - d{x)'-f^d{x)], 

il'''''{x)=u{x)-iyd{x) . 



(36) 



The correlation functions in (|35|l involve several in- 
variant functions that can be separated by the ap- 
propriate projections. Lorentz structures that are 
most useful for writing the LCSRs are usually those 
containing the maximum power of the large mo- 
mentum p+ ^ pz. We define, for the loffc current 



A+r,i'""= 



ipz){mMAr+^^Br}N+{P), 



A+T^-^={pz) {mNAt+ ii_Bl}^^N+{P) , 
K+Tl'^={pz) {m^Al+i^BDN+iP) , (37) 

where A and B depend on the Lorentz-invariants 
Q2 = -q2 ^^^ pi2 = (p _ g)2 poj. the leading- 
twist current 773 we use instead 



T. 



,3^,„,^2(pz^ ^^^ ^^ ^_^p^ 



T: 



„a 2(pz)3 



771 



{mNAl+fi^BI}j,N-iP), 



N 



T: 



'3,v^ 



2{pzf 



m 



{mNAl+4^Bl}N-{P), (38) 



N 



and similarly for 774 

4,cm_2(P^)^ 



t: 



{mNAT+ii^BT]N-{P)^ 



"-N 



T4.a=^(Pf)! {rnNAl+ i^Bl) ^,N-{P) , 



"-N 



t: 



4,v^ 



2(pz)3 



{mNAl+fi^Bl}N-iP), (39) 



'AT 



The clastic nucleon form factor contribution of in- 
terest corresponds to a pole term in the variable 
P'^. For the relevant projections wc get 
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(40) 



Note that the pseudoscalar form factor Gp does 
not contribute to the z^T^ projection because of 
the condition qz — Q. It can be extracted from the 
LCSR for other structures, or through the relation 
to the axial form factor in Eq. (|17|l which is exact 
to our accuracy. 

On the other hand, the correlation functions can 
be calculated in QCD for sufficiently large nega- 
tive -P'^ and q^ = —Q^ in terms of nucleon DAs 
using the OPE. The corresponding expressions (to 
tree-level accuracy) are collected in Appendix A. 
Matching between the two representations one ob- 
tains the light-cone sum rule that relates the nu- 
cleon form factors with nucleon DAs. The precise 
procedure was described many times in the litera- 
ture (see e.g. j^) so we omit the technical steps. 
The resulting sum rules depend on two parame- 
ters: the continuum threshold sq ~ (1.5 GeV)^ 
and Borel parameter M^ which defines the scale 
at which the matching between the two represen- 
tations is done. The dependence on the Borel pa- 
rameter is rather weak. For definiteness. in the 
plots shown below we take AP = 2 GeV^. 



C. Results: lofFe current 

In this section we present LCSR predictions for 
the nucleon form factors that are obtained using 
the loffe interpolating current 771 for the proton. 
The form factors are plotted in the range of the mo- 
mentum transfers 1 < Q^ < 10 GeV^; for smaller 
Q^ our approach is not applicable, for larger Q^ 
we expect that radiative corrections to the sum 
rules (that include in particular the usual pQCD 
contribution) will become dominant. The calcu- 
lations are done using two representative sets of 
nucleon distribution amplitudes: asymptotic DAs 
(solid curves) and including the corrections esti- 
mated using QCD sum rules (dashed curves), see 
Appendices B,C,D for the definitions. At this stage 
we do not attempt to fit the form factors by tuning 
the parameters of DAs, the difference between the 
solid and the dashed curves gives more or less the 
range of form factor values that can be obtained 



1.75 
1.5 

1.25 

1 

0.75 

0.5 

0.25 



G^/(HpGd) 



0^ 



10 



FIG. 3: LCSR prediction (loffe current) for the mag- 
netic form factor of the proton normalized to the dipole 
form factor G'j^j / {fipG d) ■ The data points: ■^: SLAC 
1994 IS; A: SLAC 1994 IH; ■: SLAC 1970 H^*; 
♦ : Bonn 1971 IH*; D : Stanford 1966 [H*; 0: JLab 
2004 0|; A: JLab 2005 0. (*: Data actually taken 
from 1531 1. 



with the DAs of "reasonable" shape. 

The prediction for the proton magnetic form 
factor normalized to the dipole form factor, 
Glf/iiipGo), where 

Gd = {I + QV0.7lGcY^f, 

is shown in Fig.O Both the Q^ dependence and the 
magnitude of the form factor is reproduced rather 
well, especially if using asymptotic DAs. 

The result for the ratio of the proton electric 
form factor to the magnetic for factor, pipG^/G^j, 
is plotted in Fig. 21 For completeness, we in- 
clude on this plot the data obtained both via 
Rosenbluth separation and the polarization trans- 
fer techniques, although the former one is most 
likely flawed. Most interestingly, this ratio ap- 
pears to be very sensitive to the shape of nucleon 
DAs. Whereas the experimental data obtained via 
Rosenbluth separation could nicely be described 
by asymptotic DAs alone, the polarization trans- 
fer data require considerable corrections. 

The LCSR predictions for the neutron are shown 
in Fig. [31 and Fig. |H1 for the magnetic and the 
electric form factors, respectively. In this case, 
again, the magnetic form factor is described rea- 
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FIG. 4: LCSR prediction (lofFe current) for the electric 
form factor of the proton normalized to the magnetic 
form factor: fipG^/C^j. The data points shown in red 
are obtained via Polarization transfer: ▲: Jefferson 
LAB 2002 Izi; ■: Jefferson LAB 2001 III; ♦: Jeffer- 
son LAB 2000 17(1 : The data points shown in blue are 
obtained via Rosenbluth separation: ■: SLAG 1994 
Ei; ♦: SLAG 1994 [H; A: SLAG 1970 [13 *; *: 
Bonn 1971 ^*; D: Stanford 1966 [H*; 0: JLab 2004 
[54; A: JLab 2005 1551 . (*: Data actually taken from 
[53 )■ (Color identification refers to the online version) 
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FIG. 5: LGSR prediction (loffe current) for the mag- 
netic form factor of the neutron normalized to the 
dipole form factor G^i/^^nGo)- The data points: ■: 
SLAG 1993 HI; ▲: Mainz 2002 [sj; ♦: Mainz 1998 



sonably well by asymptotic DAs, while the magni- 
tude and even the sign of the electric form factor 
depends on their shape. Further, the LCSR pre- 
diction for the axial form factor of the proton G^*^ 
normalized to G)y = 1.267/(1 -|- Q^)^ is shown in 
Fig.[7| Experimentally, this ratio is close to one. A 
more steep Q^-dependence of the axial form factor 
compared to the electromagnetic ones seems to be 
correctly reproduced, and also the normalization 
agrees within 50% accuracy. 

Last but not least, in Fig.|Slwe show the LCSR 
result for the form factor G^*-^ normalized to G^*^. 
As mentioned above, this form factor is forbid- 
den by T-invariance so that this ratio has to be 
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FIG. 6: LGSR prediction (loffe current) for the electric 
form factor of the neutron G^{Q^). The data points: 
Red symbols: experimental values: ♦: SLAG 1993 1561 : 
A: Jefferson Lab 2001 @; ■: Mainz 1999 0. 
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FIG. 7: LGSR prediction (loffe current) for the axial 



form factor Ga normalized to G 



(a) 



9a/{1 + Q' 



zero. In the LCSR approach the final and the ini- 
tial state nuclcons arc treated differently and the 
T-invariance is not manifest. Smallness of Gt is 
therefore an indication of how good the nucleon 
state is separated from the continuum by the simple 
duality assumption. We observe that G^'-^/G^'^ 
strongly depends on the shape of the nucleon DAs. 
It is small and negative for asymptotic DAs but be- 
comes positive if the DAs acquire large corrections. 



D. Results: Leading-twist currents 

1. Checking isospin relations 

One of the main motivations for our study is to 
find out the optimal nucleon current for the calcu- 
lation of nucleon form factors within the framework 
of LCSR. 

In Ref. [331 the current 774, Eq. |j23, was used 
which in difference to 773 , Eq. (|32|) , couples both to 
isospin / = 1/2 and / = 3/2 states. A priory, it is 
not obvious which current is better since the sum- 
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FIG. 8: LCSR prediction (lofFe current) for the tensor 
form factor GjJ normalized to G^ • 



niation over quantum numbers may improve the 
accuracy of the duahty approximation for the con- 
tinuum and also it usually suppresses poorly known 
contributions of higher-dimension(twist) operators. 
The new observation of this work is that this ar- 
gumentation can be tested by checking the isospin 
relations e.g. for the vector current: 



Ft^{Q')=Ff{Q')-F^{Q^ 



F^^iQ-')^F^iQ-')^F^\Q'), 



(41) 



cf. H14|l . These relations are fulfilled identically if 
the 773 (or rji) current is used, because in this case 
the correlation functions H35() satisfy isospin rela- 
tions by themselves. However, when using the 774 
current, the extracted form factors satisfy the re- 
lations in 1)41(1 only approximately, within the sum 
rule accuracy. In particular, their violation pro- 
vides one with a direct measure of the contamina- 
tion of the nucleon contribution by isospin / = 3/2 
states. 

The results are shown in Fig. |51 for Fi and F2 
on the two upper and two lower panels, respec- 
tively. We see that (unphysical) isospin breaking 
is relatively moderate in case that asymptotic nu- 
cleon DAs arc chosen, but it explodes if the DAs 
acquire significant corrections. The situation with 
axial form factors proves to be similar. Since the 
ultimate goal of our study is to determine nucleon 
DAs from the comparison to the data, such a be- 
havior presents a crucial disadvantage. We con- 
clude that the LCSRs with the current 774 do not 
pass the test; they are strongly contaminated with 
the isospin / = 3/2 contributions and do not allow 
for any quantitative form factor determinations. 
Hence, hereafter we drop the 774 current and only 
consider 773 (and 7/1). 



2. Nucleon form factors for rja 

Here we present LCSR results for the nucleon 
form factors, obtained using the interpolating cur- 
rent current 773. Electromagnetic form factors are 
shown in Fig.^J and weak form factors in Fig. 1111 
respectively. As above, the calculations are done 
using asymptotic DAs (solid curves) and includ- 
ing the corrections estimated using QCD sum rules 
(dashed curves), see Appendices B,C,D for details. 
The agreement with the data is in general some- 
what worse compared to the calculations using the 
loffe current and, most interestingly, the correc- 
tions to asymptotic DAs "work" in opposite direc- 
tion. We repeat that nonzero values obtained for 
the tensor form factor Gt are artifact of our ap- 
proach and can be used to quantify the error esti- 
mates. 

As a matter of principle, sum rules using all in- 
terpolating currents have to produce the same re- 
sults. In practice it has never been the case and the 
optimal choice of the interpolating current is a very 
important part of the QCD sum rule method. In 
our case it is possible that the difference between 
predictions based on 771 and 773 currents will de- 
crease when radiative corrections to the sum rules 
are included. Still, on the basis of information that 
we have now and the experience of QCD sum rule 
calculations with baryons in general, we believe 
that the loffe current 771 provides the best option 
for the construction of the LCSRs. The LCSRs 
based on the leading twist current 773 are valid and 
useful for making a consistency check since their 
structure and the relative weight of DAs of differ- 
ent twist is very different. 

To illustrate this issue we show in Fig. ^] the 
LCSR results for Ff and F|' obtained using the in- 
terpolating current r/i (upper two panels) and 773 
(lower two panels) and asymptotic nucleon DAs. 
On each plot, solid curves correspond to the sum 
of contributions of all twists and the dashed-dotted 
curves show the contribution of leading-twist DAs 
only, including the corresponding nucleon mass cor- 
rections. Notice that the Q^ dependence of the 
leading-twist and the higher-twist contributions is 
almost the same. This is to be expected, since 
higher-twist corrections to the sum rules are only 
suppressed by a power of the Borel parameter, 
not a power of Q^. On the other hand, the rel- 
ative weight of the leading-twist and the higher- 
twist terms depends strongly on the current: E.g. 
for loffe current Fi is almost entirely higher- twist, 
whereas for the 7/3 current the leading-twist con- 
tribution is dominant. Also for F2 the sum rules 
based on the 773 current are more sensitive to the 
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FIG. 9: Test of the isospin relations for the form factor Fi (upper panels) and F2 (lower panels) using the 
asymptotic form of the nucleon DAs (left panels) and including corrections estimated from QCD sum rules (right 
panels). The red (lower) and the blue (upper) solid curves show the r.h.s and the l.h.s. of the relation in 14111 . 
obtained from the LCSRs using the current 774 and normalized to the dipole form factor Gd- The blue-red thick 
dashed curves shows the LCSR result obtained with the 773 current, in which case the isospin relation is satisfied 
identically. (Color identification refers to the online version.) 



leading-twist DAs and may be more useful to ob- 
tain restrictions on the corresponding parameters. 
We expect that such differences will be moderated 
upon inclusion of the radiative corrections to the 
sum rules. However, the loffe current-based sum 
rules will most likely still provide higher accuracy 
for the form factors. 



E. A model for the nucleon distribution 
amplitudes 

The nucleon DAs provide the principal nonper- 
turbative input to the LCSRs. As we have seen, 
in many cases experimental data are in between 
the LCSR calculations that asymptotic and QCD 
sum rule-based DAs. This suggests that a good 
description of the data is possible by tuning the 
parameters of the DAs. As a demonstration, we 
present here the results obtained using a simple 
model in which the deviation from the asymptotic 
DAs is taken to be one third of that suggested by 
the QCD sum rule estimates. 

The corresponding parameters are: 



Al 



ft 



= 0.13, 
0.33, 



£U _ 
Jl ~ 



= 0.30, 
0.09, 



/I 



0.25. (42) 



where the first two refer to the leading twist-3 and 
the rest correspond to twist-4. These values are 



not unreasonable, since QCD sum rules are known 
to overestimate the matrix elements of higher con- 
formal spin operators, and we just made the sim- 
plest assumption that all sum rule results have to 
be rescaled by the same factor. Our leading-twist 
parameter Vi is very close to the phcnomenologi- 
cal Bolz-KroU model J93; Ai is somewhat bigger 
but the dependence of the leading-order sum rules 
on this parameter is weak. To this accuracy, the 
sum rules also do not depend on the parameters A2 
and /^; this dependence is present, however, in the 
transition form factors like ■y*N -^ A. 

The calculations using this model are shown by 
solid curves for the electromagnetic form factors in 
Fig. El and weak form factors in Fig. E) respec- 
tively. In addition, in Fig. 1151 we plot the corre- 
sponding y/Q^F2 /Ff ratio. On the same plots we 
show by the dotted curves the effect of the variation 
of the couplings ratio fw/^i (see Eqs. (1^ within 
a conservative 30% error range. This ratio deter- 
mines the overall normalization of the leading-twist 
DAs compared to higher twist so that the sensitiv- 
ity to /jv/Ai is a good indication of the relative size 
of the leading- twist contributions to the LCSRs. 

One sees that the experimental data on the elec- 
tromagnetic form factors are reproduced very well, 
and, most welcome, the unphysical tensor form fac- 
tor Gt becomes consistent with zero. Also for the 
axial form factor there is a good agreement, both 
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FIG. 10: LCSR results for the electromagnetic form factors of the nucleon, obtained using the leading-twist-3 
interpolating current rj^. Identification of the curves and the data points is the same as in Figs. 13101 




FIG. 11: LCSR results for the axial form factor of the proton GJ normahzed to G^ = Qa/O- + Q^)^ (left 
panel) and tensor form factor Gj- normalized Ga (right panel), obtained using the leading-twist-3 interpolating 
current 773. 



in shape and normalization. 

Last but not least, we can use the same set of 
DAs to calculate the ^*N -^ A transition form fac- 
tors within the LCSR approach, following ReL |4g. 
The results are shown in Fig. ^] In this case we 
also get a much better agreement with the exper- 
imental data on the electric form factor compared 
to the calculations that use asymptotic or sum rule- 
based DAs. 

We should warn that the model in Eq. H42|) is not 
based on any systematic attempt to fit the data and 
in fact we believe that such any fitting would be 



premature before the radiative corrections to the 
LCSR are calculated. In addition, one has to take 
into account the scale dependence of the parame- 
ters of the DAs and study in more detail the de- 
pendence of the sum rules on the Borel parameter. 
Still, the very possibility to describe many different 
form factors using the same set of DAs is nontrivial 
and indicates the sclfconsistency of our approach. 
The true parameters of the DAs are probably not 
far from the numbers quoted in Eq. (|42|l . although 
at this stage we cannot give any error estimates. 
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FIG. 12: LCSR results (solid curves) for the electromagnetic form factors of the nucleon, F[ (left) and F|' (right) 
obtained using the interpolating current rji (upper two panels) and 773 (lower two panels) and asymptotic nucleon 
DAs. On each plot, the dashed-dotted curves show the contribution of leading-twist DAs only, including the 
corresponding nucleon mass corrections. 
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FIG. 13: LCSR results (solid curves) for the electromagnetic form factors of the nucleon, obtained using the model 
of the nucleon DAs I42I I and the loffe current rji . The dotted curves show the effect of the variation of the ratio 
/at/Ai by 30%. Identification of the data points is the same as in Figs.|2HSl 



IV. CONCLUSIONS 



Following the proposal in Ref. |35J we have made 
a comprehensive study of leading order light-cone 
sum rules for the electromagnetic and weak nucleon 



form factors. We presented detailed results that are 
obtained using different interpolating currents for 
the nucleon and argue that the loffe current ap- 
pears to be the optimal. We make an update of 
the QCD sum rule estimates of the shape parame- 
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FIG. 14: LCSR results (solid curves) for the axial form factor of the proton G5 normalized to G^ = <7A/(l-fQ^)'^ 
(left panel) and tensor form factor Gy normalized to G^ (right panel), obtained using obtained using the model 
of the nucleon DAs l|42|l and the lofFe current ?7i . The dotted curves show the effect of the variation of the ratio 
/at/Ai by 30%. 
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FIG. 15: LCSR results (solid curves) for the ratio 
\/Q^F2 /Fl obtained using the model of the nucleon 
DAs 1421 1 and the loffe current rji. Data points: Red 
symbols: experimental values obtained via Polarization 
transfer: A,^: M. Jones (private communication); 
Blue symbols: experimental values obtained via Rosen- 
bluth separation: ■: SLAG 1994 [H; ♦: SLAG 1994 
[431 . (Golor identification refers to the online version) 



ters of higher-twist DAs and also present a simple 
model that seems to reproduce main features of the 
experimental data remarkably well. In addition we 
included a complete summary of higher-twist DAs 
and some new expressions for the OPE of three- 
quark operators which extend the results given in 

Our main conclusion in this work is that the 
LCSR approach to baryon form factors seems to 
be sufficiently accurate to allow one to get a quan- 
titative description of hard exclusive reactions with 
baryons. From the theory point of view, this tech- 
nique is attractive because in LCSRs "soft" contri- 
butions to the form factors are calculated in terms 
of the same DAs that enter the pQCD calculation 
and there is no double counting. The asymptotic 
pQCD limit, in fact, formally corresponds to a part 



of the two-loop a^ radiative correction to the sum 
rules. Thus, the LCSRs provide one with the most 
direct relation of the hadron form factors and dis- 
tribution amplitudes that is available at present, 
with no other nonpcrturbative parameters. 

We remind that the sum rules considered in this 
work are tree-level. Further progress in this di- 
rection requires the calculation of the radiative 
one-loop correction to contributions of twist-3 and 
twist-4 operators. It would also be very valuable 
to have lattice evaluations for at least some of the 
parameters that enter the DAs, most importantly 
the /at/Ai ratio. 
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APPENDICES 

APPENDIX A: SUMMARY OF 
CORRELATION FUNCTIONS 

In this Appendix we present the tree-level results 
for the correlation functions defined in Eq. (|37|) . 
(|38|) . The correlation functions are expressed in 
terms of nucleon DAs that are summarized in Ap- 
pendix B below. We use the following notations: 
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FIG. 16: 7 * A'^ — > A transition form factors in the LCSR approach |4y]. The blue solid curve corresponds to the 
calculation using the model in Eq. I42II . The red dash-dotted and dashed curves are obtained using the asymptotic 
and the QCD sum rule motivated DAs, respectively. For the references to the experimental data see |4y]. (Color 
identification refers to the online version) 



rl — x^ 

I nil 

X / dx\F{x\^\ — x\— X'^^x-^^ (A.l) 
Jo 

and 

_ rxi /-l-aJa 

F{x2)^ dx'2 j dxiF{xi,x'2,l~ xi- x'2) , 

F{X2) ~ / dx'2 / dX2 

r^-^-2 „ 

X / dxiF{xiTX2T^ ~ xi — X ), (A. 2) 

JQ 

where F = A,V,T \s a. generic nuclcon DA that de- 
pends on the three valence quark momentum frac- 
tions, and also shorthand notations for the combi- 
nations of the DAs: 

F123 ^Vi~V2-V^, 

1^1345 = -2^1 + V^ + Vi + 2^5 , 
F12345 - 2^1 - F2 - 1^3 - 14 - ^5 , 



V^ 



123456 



-Vi + V2 + Va + Vi + V5~V(i, (A.3) 



^34 = A3 - ^4 , 

^123 = -Ai + A2 - ^3 , 

^1345 = -2Ai - ^3 - A4 + 2A5 , 

^12345 = 2^1 - ^2 + A3 + ^4 - ^5 , 

Ai23456 = Ai - A2 + A3 + A4 - As + As (A.4) 

and also 

^137 ~ T1 — T3 — T-j , 
13478 = -^-Ll — -1^3 — -L4 — -L7 — -L8 , 

^134678 = Ti-n-Ti + Te-Tj-Ts. (A.5) 

In addition, in this Appendix we use 

q3 = q- X3P , q2 = q- X2P ■ 

In this notation we obtain, for the correlation 
functions involving the loffe current 771: 
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AT= 2e, /' dx, (^-t^yi23 + 5 r dx,V^{x,) + ^^^43 



3 1 4 

I 93 



4 Jo 



2e„ ! dx^l^ I ' dxi [-2Vi + 3^3 + A3] i^r) - 
Jq I 92 ^0 






r-qi 



Vv. 



—^ •'123 

9| 



^2^2 r 



Q + 92 ;r ^^"^N 

4 ^123 4 

92 92 



V1345 - 2^43 + A 



34 



2x2m% 



Vi 



23456 



Jo 1 93 Jo 






9| 



4 



T423 - "^^4 



2eti / dx2 
Jo 



1 



t2 r 

N 



dx,[V^+A^]ix.)+■^ V.^'^'+At 



X2mjf "■ 
9| 



92 Jo 92 

1^345 + yi23 + A123 — 2V43 + ^34 



,M(u) ^ .Miu) 



(A.6) 



At 



2 I „2 



Q^+93' 



93 



X3_ 

9I Jo 



-c / rf^.^^:^_j«^^^^^:::| / rf^^v^3(^^) + r:3:^^^3 



9 2 

XoTni 



93 



rfa-2 



2^2 



2 s — / dxi[2Ai+ 3 A3 + 1/3] {x^ ) + ^^4^ Ai 
I 92 Jo 93 



2a;2"% ,A/(«) , Q - 92 :? 

4 ^123 



93 



- 1^123 



9| 



9| 



^1345 ~ 2^34 + 1/43 



2x2m-j^ 



A12 



4 ^^123456 
93 



6^ 



dx3 



93 Jo 93 



4 



^123 — ^3 



f 1 f^^ 

dx_ 
I 92 Jo 



2 ^ ^ / dxi [Fi + Ai] (x,) + ^ Vi'^ ^"^ + Af 



,M(«) , -Af(«) 



X2ml, 

4 



^1345 ~ ^123 ^ ^123 " 2A34 + V43 



(A.7) 



For the correlation functions involving the leading-twist current 7^3 we get 



AT=--eu / d.T2 <— dxi 
-J Jo I 92 Jo 






93 



^2"'-N 

~4~ 



^12345 + 2?! 



12345 ~r ^J 13478 



ZX2TTlj^ 



3 Jo I 92 



^123456 ~ 2Ti34678 

V123 + 2Ti37 



- -^ed / Axz\^X2 ^ xz,F ^ F)^ 






1^123456 — 2Ti34678 



^ed / da;3|(x2 ^ 2;3),F^ i^l 



(A.; 



Ar'=\ fdx2 (^ /''d^i Ai + 2Ti] (x,) + ^ 
-J Jo I 92 Jo 92 



^j\/(«) _^ 27-A^(") 



''2"''N 



j4i2345 + 2ri3478 






A 



123456 + 2Ti34678 



3 Jo ^n 9i 7o 



dx3<i9 / dxi[Vi{x,)~2Ti{x,)] 
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X3m% 

' 4 



\^M{d) _ ^^M(d) 



4 



ZiXnin j<j 



4 



Bl 



3 Jo 



1 



dxo 



dx- 



X2m% 

4 



I 4 



A 



123 



2Ti 



137 



V12345 ^ 2ri3478 

^123456 + 2Ti34678 



1^123456 + 2Ti34678 



V. 



123 



2Ti 



137 



Zx2Tnj^ 

ZX^TTlj^ 



4 



3 JO 

and, finally, for the correlation functions involving r/4: 
»i 



F 



123456 



2T 



134678 



da:2 







X2 



Q2 Jo 



dxiVi{Xi) H 74— Vj 74— V^12345 



?l 



9l 






V 



123456 



(A.9) 



2erf/ dx3{{x2^X3,F->F)'j, 



Br= 4e„ / d.T2 



a:2W^ 
9l 



2^4 



^^123 + 



2xim 



2'"- AT 



9| 



y 



123456 



2ed da;3{(a;2-*a;3),F^F}, (A.IO) 



Ar 



2 /dxa 



dx. 



a;2 



92 Jo 



X3 



dxiAi{xi) 



X2m% A/(u) xlm% ^ 

4 — -^1 4 — ^12345 

4 4 



3^4 



2.TSm 



2"'-Na 



„2^2 



. TA ^ \ , xsm^ M{d) Xy,m.^ ~ 

3\— I dxiVi(x^) + ^Vi ^Vi 

93 Jo % 



4 



2345 



9| 



123456 



4 



V 



123456 



s: 



K,2^4 



^2 / dx-2 "^ 4^^123 



4 



2'"JV ^ 

— S ^^123456 

92 



da:;3 



4 



Vi 



123 



4 



V 



123456 



(A.ll) 



AV'^-Z I d:r2 



L 92 Jo 



''2'"jV 



Br'^2 / dx2 



X2 r . , , s X2mtj 

4 / ^2:1 [l/i - ^1 - 2Ti] (xO + ^^ 

92 



2 r 



M(«) .M(«) r,-r^-'"(") 



1^2345 + ^12345 + 2T13478 
^123 + ^123 + 2Ti37 



X2m\ 



4 



^123456 + ^123456 + 2Ti34678 
^123456 + ^123456 + 27x34678 



,(A.12) 



AT'=2 /dx2 (^ Cdx^ [Fi - Ai + 2Ti] (x,) 
Jo I 92 Jo 



2 2 

4 



BV'=2 /"'dxs/^^ 



V 



1^2345 + ^12345 ^ 2T13478 
^123 + ^123 ~ 2ri37 



4 



X2m%^ 

^ 4 

4 



yAf («) _ _4^A/(«) ^ 2T^'^'^'> 



V 



+ ^1 



123456 "T ^123456 



2r 



134678 



V 



123456 + ^123456 



2T 



134678 



,(A.13) 



In all expressions the functions with a "tilde" and a "hat" have X3 and X2 as an argument, respectively, cf. 
(jA.ip . 1A.2|) . Also, in the terms involving two integrations over the momentum fractions, the remaining 
momentum fraction is replaced by using xi + a;2 + 0:3 = 1. The results in (jA.lOp agree with the cor- 
responding expressions in |35( up to two misprints: a factor two in the V'^^'^^-term and the sign of the 
^123456 contribution; the other expressions are new. 

The answers for the neutral vector current JZ''^^ are easily obtained from the corresponding expressions 
for the electromagnetic current by a substitution e„ -^ 1/2 and e^ -^ ^1/2- Also, since the currents 771 
and 773 are pure isospin / = 1/2, the correlation functions involving the flavor-changing charged currents 



18 



j^''^'^ and j^''^'^ are given in terms of the corresponding correlation functions involving the flavor-conserving 
currents by exact isospin relations 



A^ 



2a: 



b: 



2SJ 



Al 



2Af 



Bf 



2Bf 



1,3. (A.14) 



These relations are not manifest because they involve isospin relations between different nucleon DAs and 
provide a nontrivial check of the calculation. 

The Borel transformation and the continuum subtraction are performed by using the following substi- 
tution rules: 



dx 



dx 



dx 



g{x) 



{q-xPf 
g{x) 



dx g{x) 
dx g{x) 



^ dx 



g(x) exp 



xQ^ 



xra 



N 



{q-xPy Jo a;2 (s - P'2)2 

g{x) /■ dx g{x) 

[q ~ xPf "" X ^ (s - P'2) 



1 /■' dx 
T72 / — e(a;)exp 



M^ J.o ^ 



'/2\3 



1 dx 



g{xo)e 



-so/M^ 



2AP ./.„ X 



1 



2xo{Q'^+xlm%)M^ 2Q^ + xlm 



N 



xAP 


AP J 


1 ' 


( xQ^ 


xrnjf 
AP 


h 


V xAP 


)exp - 
V 


xQ' 
xAP 


AP 


■ d 


q{xo) 





g{xo)e 



-so/M^ 



2 2 
XqUIj^ 



dxo xo ((52 + xlmjj) 



,-so/M^ 



(A.15) 



where M is the Borel parameter, s = ^-^Q^ + (1 — x)m'j^ and xq is the solution of the corresponding 
quadratic equation for s = sq: 



Xo = 



y^(Q2 + .so-m^)2+4m^02 _ (q2 ^ ,,^ _ ^2^^ 



I {2ml) ■ 



(A.16) 



The contributions ~ q~'^o/m jj^ gq^ ljA.15p correspond to the "surface terms" arising from successive 
partial integrations to reduce the power in the denominators {q — xPY^ = (s — P' )2^(— a;)2^ with 
A^ > 1 to the usual dispersion representation with the denominator ^ {s ^ P' ). Without continuum 
subtraction, i.e. in the limit sq ^ 00 these terms vanish. 

In addition, in the hadronic representation for the same correlation functions one has to make the 
substitution 



1 



JM^ 



m 



N 



pi2 



(A.17) 



As an example, we present here the final sum rules for Ff^ F2 and G^ obtained using the lofFe current: 



icni //o2\ 



Fr(Q 



2Ai 



lein //-j2 N 



FriQ 



GT{Q')=i^^ 



dx 



dx 



dx 



Q2{x) , gl{x) 
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xHP 

q\{x) ^ g\{x) 
X x'^A'P 
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X x2M2 



exp 



exp 



exp 



xQ'^ xmlj' 
xAP AP , 


\ (?|(a;o)e-(""-™«)/*^' 


j ' Q^ + xlml 


xg2 a;m^\ 
xAP AP ) 


e|(x-o)e-(^«^™«)/*^'' 


(32 ^ xlml 


xQ"^ xm^' 


\ e^(xo)e-(««-"«)/*^' 


xNP ' A/P , 


1 ' Q^ + xlml 



(A.18) 



with 



(?^(x)== 2ed{ 1^123 +x dxiV3.{x^) +\+2eulx dxi [-2^1 + 8^3 + ^3] ix^) - V123 + ^123 
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glix)^ 2eA QVi23 + x^m%Vi3 >+2eu{ Q^ F123 + ^123 - x^m% 



^1345 — 2V43 + ^34 



2a:m^ ( Vf '"^ + ^123456 



=-2ed i / dxiViixi) I + 2e„ | / dxi [Vi + Ai] (xi) 



+ X i>i345 + 1^123 + ^123 " 2^43 + I34 

^2(2;)=! V123 +x dxiVz{xi) \ + \x dxi [2Ai + 3A3 + V3] (xi) - A123 + V123 



Ql{x) = i Q^Vi23 + a;^TO^^43 \ + {Q [ ^123 + ^^123 + x'm 



'■N 



^1345 ~ 2^34 + V43 



2xm% A^ ^"■' -^123456 



(A.19) 



Form factors of the neutron are obtained by the substitution e^ <-> ea- 

APPENDIX B: NUCLEON DISTRIBUTION AMPLITUDES 

In the foUowing we give a summary of the three-quark distribution amphtudes from twist-3 to twist-6 
as obtained in |83l | . The DAs are defined by the matrix element of the renormahzed three-quark operator 
at hght-hke separations 

4 (0| £^^''=<(aiz)«^^(a2z)d^(a3^) 1^) = 

=S'imjvCa/j (75A^+)^ + S2mNCo,ii (tsA^")^ + PiruN (tsC)^^ N+ + P2mw (tsC)^^ A^" 



7 2 



V4 , „, /I __N ,, m^ 



+^m^ (7^C)„, (t-'tsA^-), + ^5^ {^C)^p h,N+)^ + ^V, (/C)„, {^,N-) 



2pz 



2pz 



A3 



+Ai {jh^C)^p N+ + A2 (^l^C)^p N~ + -^mi, (7±75C)„^ {i^N+) 



A 



+^rnN {l±l,C)^p h^N-) -f A5-f {^l,C)^p N+ + -^A^ {^l,C)^p N. 



2pz 



2pz 



TUN 



+Ti {icT^,C)^p (7^75A+) + T2 {io^,C)^p (7^75A-) + T3^ {^^p.C)^^ {l^N+)^ 

pz 



TUN 



+T4^ (^a.,C)„, (75A^-), + Ta^ (z.^.C)„, (7-^75A^-^), + '^J, (^.^.C)„, {^^^,N-), 



n 



+m^y (^±±'C)„^ (cT^^'75A^+) +m^^ (^±±'C)„^, (a^^'75A^-) 



(B.l) 



where for brevity we do not show the Wilson lines that make this operator gauge-invariant; a, /3,7 are 
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Dirac indices and C is the charge-conjugation ma- 
trix 113. Each DA F = V^,A,,Ti,Si,Pi can be 
represented as 

F{aj,Px)= /l?a;e-'^^^-^-"-i^(a;,), (B.2) 

where the functions F[xi) depend on the dinien- 
sionless variables Xi,Q < Xi < I, ^ ^ Xi = 1 which 
correspond to the longitudinal niomcntum frac- 
tions carried by the quarks inside the nuclcon. The 
integration measure is defined as 



h^L 



da;ida;2da;3 8{xi + X2+ X3, — 1) . (B.3) 



Distribution amplitudes can be expanded in con- 
tributions of operators with a given conformal spin 
[83| . This is convenient since operators with differ- 
ent spin do not mix under rcnormalization in one 
loop. More importantly, only the operators with 
the same spin can be related by equations of mo- 
tion so that the truncation of the conformal spin 
expansion at a certain order produces a sclfconsis- 
tent approximation. In Ref. [83| the contributions 
of the leading and the next-to- leading order in con- 
formal expansion have been taken into account. 

To this accuracy one obtains twist-3 DAs: 

Vi{x,, ^i)=l2QxiX2X:i [(/.!](//) + 0+(//)(l - 8x3)] , 

Ai{Xi, ^i) = l2QxiX2X3{x2 - Xi)(l)^{ll) , 

Ti[xi,y)=l2QxiX2X3, 4>%{y) 



(B.4) 



Twist-4 distribution amplitudes: 

V2{X„ ^i)^2AxlX2 [:fM + <^t(M)(l - 5x3)] , 

A2{xi, iJ.)=2AxiX2{x2 -a:i)^4 (^), 

T2ix„ fi)=24xiX2 [e2(M) +e|(/^)(l - 5x3)] , 

V^ix,, fi)=12x3[i''i{li)il - X3) 

+ ip^{n){l ~ X3~ 10a;ia;2) 

+ip^{n){xj + x\ - xj,(\ - x-i))\ , 

A3{x,,ij)=12x3{x2 - xi) [{i^t+i^t) (a*) 

+ ij^ {fi){l - 2x3)] , 
T3(x„m)=6x-3 [{<Pl + V!^ + e2)(M)(l - X3) 

+ i^t +^++ (t)ifi)il -X3- 10x1x2) 

-t-((/)4 - ■(/'I + 0(A«)(2;i +xl- X3(l - X3))] , 



Tr{x,,^i)^6x3 [{cf>1 + ^l - ^4°)(m)(1 - 2:3) 

+ {(t>t + V^l - C|)(m)(1 - 2^3 - 10X1X2) 

+{4>A - V'l - S,'^){^J){x\ + xj - X3(i - X3))] , 

5i(x,, /i)=6x3(x2 - xi) [{<j)l + i^l + ^2 + (t>t + i^t 

Ht)i^^) + (<^r - ^r + 0(/^)(i - 2x3)] , 

Pl(x,, aO=6x3(xi - X2) [((/.^ + i'°4 - ^4 + <l>t + ^4 

-C){^^) + ('^r - ^I - 0(a*)(i - 2x3)] . 

(B.5) 

Note that T3 and Ty differ only in the sign of the 
^-contributions, while Pi and ^i differ only in the 
sign of the </>- and ■i/'-contributions. 
Twist-5 distribution amplitudes: 

V4ix,,^l)=3[lpU^^)il~x3) 

+ ^+{fi){l-X3-2{xl+xl)) 

+ i^sif^) (2X1X2 - X3(l - X3))] , 
A4(X,;, Ai)=3(x2 - Xi) [-^bUl^) + ^^(/X)(l - 2x3) 

Ux„^,)=^ [(0° + Vs" + e5°)(M)(i - 2:3) 
+ i^t + i't + U) (a^)(i - ^3 - 2(x? + xD) 

+ {^b - i-^5 + ^5) (m) (2xiX2 - X3(l - X3))] , 

r8(x„M)=^[(0g + V^,?-e5°)(M)(i-^3) 
+ i^t + i't - 4+) (a^)(i - ^3 - 2(x? + xD) 

+ i^l^b - ^5 - ^5) (m) (2xiX2 - X3(l - X3))] , 
y5(x„M)=6x3 [</'g(M) +0^(/i)(l - 2x3)] , 
A(Xi,M)=6x3(x2 -Xi)(p^{fl), 
Tslx,;, m)=6x3 [^lifl) + e5+(M)(l - 2x3)] , 

52(X„ m) = ^(x2 - Xi) [- (0° + ^5° + eg) (/^) 
+ (0^+^'5++e5+)(/^)(l-2x3) 

+ i'l^b -^5 +^5) (m)2;3] , 

P2(X„ m) = ^(x2 - Xi) [- (-0^ - V^g + ^5") (m) 



+ M^-^5++e5+)(M)(l-2x3) 



(B.6) 
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Note that T4 and Tg differ only in the sign of the 
^-contributions, while P2 and S2 differ only in the 
sign of the (p- and i/j-contributions. Note that the 
results for S2 and P2 quoted in [83j contain mis- 
prints. 

Finally, twist-6 distribution amplitudes: 

Veix,,fi)^2 [0g(/i) + 4(/.)(l - 3x3)] , 

Ae{xi, fi)=2{x2 - xi)(I)q , (B.7) 



T6(x„M)=2[0^(Ai) - 2 (4 - ^e) (l-Sxg) 



In all cases fj. is the renormalization scale. 

The coefficients in the above expansions can 
be expressed in terms of eight non-perturbative 
parameters Jn, Ai, A2, /f, ff, fi, A^A^f- which are 
defined in Appendix D, see also [23 j section 3.2. 
The corresponding relations read, for the leading 
conformal spin: 



^t^^^fN, 0^ = 01? = ^(./n + Ai), 



^4" = e" = 7jA2 , 7^4" = V^° = ± (/^ - Ai) .(B.8 



1 
2 
1 
2 



For the next-to-leading spin, for twist-3: 

03 = f /iv A?, 0+ == ^/^ (1 - 3V,^), (B.9) 
for twist-4: 
<A| = J [./w(3 - lOVi'^) + Ai(3 - 10/f)] , 

04-=-^ [/a^(1 - 2A^) - Ai(l - 2ft - 4/r)] , 
^4+ — J [/iv(2 + 5A5'-5Fi'^) - Ai(2-5A^-5/r)] , 
^4" = 5 [/w(2 - ^5" - 3Fi'^) - Ai(2 - 7/f + A")] , 



C = Y^A2(4-15/2''), ^4^ = -A2(4-15/2''),(B.10) 
for twist- 5: 

<t>^=-l[fN{l-2A-)-X^{ft-fr)], 
^5^~l [fN{5 + 2A-,-2V,^) - Ai(l-2/{'-2/n] , 
i^5^ I [/a^(2 - ^1 - 31^1'^) + Ai(/f - /n] , 
^5+= Aa2(2-9/|), e,^=-^A2/2^ (B.ll) 



and for twist-6: 

K^l [fNil - 4n') - Ai(l - 2ft)] , (B.12) 

06 = ^ [/^(l + 4^?) + Ai(l - 4/1^ - 2/r)] . 

The normalization of all DAs is determined by 
three dimensionful parameters fN, ^1,^2 that are 
well known from the QCD sum rule literature and 
correspond to nucleon couplings to the existing 
three different three-quark local operators with the 
correct spin and isospin, see [8J|. The numerical 
values (at the scale /x = 1 GeV) are |83ll88|: 

fN={5.0±0.5) ■ IQ-^GcY^ 
Ai=-(2.7±0.9)-10"2QgV^ 

A2 = (5.4 ± 1.9) • 10-^GeV^ (B.13) 

see also Appendix D. The remaining five param- 
eters determine the shape of the DAs (deviation 
from the asymptotic form) and their values are 
much more controversial. Asymptotic DAs corre- 
spond to the choice 



il - -^Q , il - ^Q 



f^- 
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(B.14) 



The leading-twist-3 parameters Vi and A" were 
calculated usin g QC D sum rules in Refs |88ll9lll92l | 
with the result |93l 



A'l = 0.38 ±0.15, 
Vf = 0.23 ± 0.03 , 



(B.15) 



while the remaining twist-4 shape parameters 
f±j,fijf2 were estimated by the same method in 
|83|. In this work we reconsider the correspond- 
ing sum rules (see Appendix D) and obtain a new 
estimate 



ft = 0.40 ± 0.05 , 
f^ = 0.22 ± 0.05 , 
/i" = 0.07 ± 0.05 . 



(B.16) 



The set of nucleon DAs obtained using the param- 

is sometimes referred to 



eters in (|B.15|I and l|B.l( 
as the Chernyak-Zhitnitsky-like model of the DAs. 
Alternatively, there exists a phenomenological 
model for the leading-twist DA [93 which was ob- 
tained by modelling the soft contribution by a con- 
volution of light-cone wave functions. This model 
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is much closer to the asymptotic DA com.pared to nique are not available, 
the CZ-modcl and corresponds to the choice [93j 

. „ , In the main text we suggest one more model that 

^ ' ' allows one to obtain good agreement with the data 

yd — 13/42 (B.17) oi^ the nuclcon form factors within the LCSR ap- 

proach. The corresponding parameters are given 
Estimates of the higher- twist DAs in the same tech- in Eq. (|42() . 

I 



APPENDIX C: OPERATOR PRODUCT EXPANSION OF THREE-QUARK CURRENTS 

In this Appendix we present the tree-level expansion of the nucleon matrix element of the three-quark 
operator in terms of nucleon DAs to the twist-5 accuracy. The general Lorentz decomposition reads p^ 

4 (0| e'^''u'^{aixyi,{a2x)d';{a3x) \P) 

+ (Vi + ^!^Vi^^) {rC)^p (75A^), + V^ruM (/'C)„^ (^ts^V), + Vgm^ (7^C)„^ (t^Ts^V), 
+V^ml {^C)^p (75iV)^ + V,m% h^C)^^ {in^^'^x,-i,N)^ + Vgrn^ {^C)^^ i^l^N)^ 

+Aiml (^75C)„^j N^ + A5m% (7m75C)„/3 (iaf'^x^N)^ + Aem^ {^-/sC)^^ (^N)^ 
+ (t, + ^!^r/^^ (P^*^,,.C)^^ {ri5N\ + %mt, (x^^P^za,,,C)^^ {^,N\ 

+%m% {x^^P^KJ^.C)^;, i^l^N)^ + Tjml (a^.C)„^ [a^" ^-f^N)^ + %m], (xV^.C)„^ (a^^x,75A^)^ . 

(C.l) 

where it is assumed that the "calligrafic" functions depend on x"^ at most logarithmically. Leaving aside 
the terms in x^, V*^,^*^ and T^^ , the rest of the functions can be expressed, to the tree-level accuracy, 
in terms of the nucleon DAs at the renormalization scale /i^ ^ \i/x'^\: 

Si=Si, 2PxS2 = Si-S2, Vi=Pi, 2PxV2 = P2-Pi, Vi = Fi, 2PxV2 ^Vi - V2 -V3 , 
2V3 = F3, iPxV4^-2Vi + V3 + V4 + 2V5, 4PxV5 = 14-^3, 
A{PxfVe = -Vi+V2 + V3 + V4 + V5-Ve, Ai = Ai , 2PxA2 = -Ai + A2 - A3 , 
2A3 = A3 , APxAa = -2Ai - A3 - A4 + 2A5 , APxAs = A3 - A4 , 

4 {Pxf A = ^1 - ^2 + A3 + A4 - As + As , Ti = Ti , 2PxT2 = Ti -K T2 - 2r3 , 2r3 = T7 , 
2Pxr4 = Ti-T2- 2T7 , 2PXT5 = -Ti +T5 + 2TS, 4 {Pxf Te = 2T2 - 2T3 - 2T4 + 2T5 + 2T-, + 2T^ , 
APxTr^Tj^Ts, A{PxfTs = -Ti+T2+T5-Te + 2T- + 2Ts. (C.2) 

I 

In the following we present the calculation of the the three-quark operator in Eq. (jC.ip in a simpli- 
0{x^) corrections to the light-cone expansion of ficd situation, where positions of two of the three 



quarks coincide. This approximation is sufficient 
for the derivation of LCSRs to the tree-level accu- 
racy. The strategy is based on the approach devel- 
oped in [9J,|93,|93- We present a slightly more ex- 
panded derivation than it was given in [35|, where 
only V*^ was obtained. Note that considering the 
vector and axial-vector Lorentz projections is suffi- 



J 
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cient since the tensor ones can be determined with 
the help of isospin relations [83j . 

Consider first the case where positions of the two 
M-quarks coincide. This situation occurs when the 
d-quark interacts with the (electromagnetic) probe, 
hence we refer to it as the d-quark contribution: 



x" (0|e*J'= [u'C-io,u'\ (0)d;(x) |P)=-x° 



Vi + ^Vf^f"^)) P„ (75A^), + V^m^P^ (^75iV), 



+ Vam^v (7a 75^)-, + ViVilfXa {15N) + Vem^Xa (;fc75iV)^ 



„a /f^i ijk \„i 



[v}C^o.l^u=]iS))dt^{x)\P)=-x^ 



(a, + ^!j^^f ('^)') P„ (^57V)^ + A^ruMPo. {^1,N)^ 



■ AzniN irfalbN) + AimlfXa {15N) + Aem%Xa (MsN) 



(C.3) 



r 



We remind that Vi,Ai start at leading twist- 
3, and hence V^ ,Ai are of twist-5. Strictly 
speaking, since we are not taking into account 
twist-6 contributions induced by 0{x^) corrections 
to V2,A2 and V3,A3, in order to be consistent we 
have to discard the contribution of Vg and Ae alto- 
gether. This contribution to the sum rules appears 
to be numerically negligible, however. 

For definitencsS; consider the vector projection. 
The meaning of the separation between Vi and 
V]^ is most easily understood upon the short 
distance expansion x^ -^ 0. In this way, the nonlo- 



cal "light-ray" operator in the l.h.s. of Eq. ljC.3p is 
Taylor-expanded in a series of local operators with 
three quark fields and the increasing number of (co- 
variant) derivatives acting on the d-quark. The 
separation of the leading twist part of each local 
operator corresponds to the symmetrisation over 
all Lorentz indices and the subtraction of traces. 
Without loss of generality, we can consider the ma- 
trix element contracted with an additional factor 
Xa, see Eq. (|C.3ll . so that the symmetrisation is 
achieved. To subtract the traces, we formally write 



»c^(a;)lit = Y1 



n=0 



I i: 



■ S/JiMj 



.J 



d"' ...d^'-diO), (C.4) 



where 'It' stands for the leading- twist part. Ob- 
serving that ^j^ = /p dt i" the subtracted contri- 
butions 0{x^) can be reassembled in the form of a 
non-local string operator: 



^2 ri 



dt- 



4 ./o dxadx 



-{Q\e^^^[u^C^u^]{Q)d^^{tx)\P) 



|e^^'=[u^CA'](0)4(x)|P) = 

;o|[£'^^'^[^^cA^](o)4(x)],_jp) 



Alternativ ely, the same result can be obtained by 

observing [9J| that the leading-twist light-ray oper- 

(C.5) ator has to satisfy the homogeneous Laplace equa- 
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tion 

^^ (01 [e'^^n-^^ A^l (OKW],t 1^) = 0- 

(C.6) 
Using QCD equations of motion the third hne in 
Eq. HC.5|) can be simpHfied to 



92 



—e'^*' [u'Cjtu^] (0)4 (te) 



dxadx'^ 

=2te'^^ [u'C-i'^u^] {Q)D^d':^{tx)+g\uons 

=2tdo,e'^^ [u''C7"m^] (0)4(ia;) + gluons, (C.7) 

where da is a derivative with respect to the over- 
all translation [9J|; for the matrix element we can 
make the substitution da — > ~iPa- Inserting this 
result in Eq. HC.5(I we finally obtain 

(0|e'^'^[w'CA^](0)4(a;)|P) = 

= (G| [e^-'"'^ [u^C^u^] (0)4(.t)]^^ \P) + ^(-z2P„) 

dtt (0| e'^^ [u'C-f°'u^] (O)d^(te) |P) 



giuons . 



(C.8) 



Notice that the r.h.s. only involves (up to correc- 
tions with additional giuons) the already known 
distribution amplitudes. This equation therefore 
allows us to determine V*^'-''-' — which appears on 
the l.h.s. of Eq. (|(L8|I — up to gluonic corrections. 
Consider the first term on the r.h.s. of Eq. I|('.8|) . 
We can write 

(0|[£^^'^-[«»C^«^](0)4(a;)]jP) = 



= - lvx[ 



iP-XX3 



(Pa;)];^Vi(75iV)^ 



from the expression given in |9; 

\ -iP-xxsl 



„2 2 2 pi 
iP-xX3 . X "i-N-^Z 



dte-'-f"-^^^*. 



by taking the derivative with respect to X3 and with 
respect to P^. One gets 



[e-^--^(P.x)],^ = {Px) 



-iP-xx^ 



2 9 2 

x^m%xi 



^^^-tP-xx^t 



[e-^--^]^^=^ 



/o 

-iPxx^ 



(C.IO) 



,-^222 

+ -a; inj^x^ 



dire 



2 —iP-xx^t 







2 P^3^ 



(C.ll) 



"2 / dtte-'^'-''=^K 
'0 

Note that we have corrected a misprint in [23 in 
the last term of the second equation, where a factor 
1/2 arises instead of 1/4. 

The corresponding contribution to Vj^ is pro- 
portional to the nucleon mass squared and involves 
the leading twist distribution amplitude, being an 
exact analogue of the Nachtmann power suppressed 
correction in deep inelastic scattering. The second 
contribution on the r.h.s. in Eq. (|C.8|) is special for 
the exclusive kinematics since it involves a deriva- 
tive over the total translation that vanishes for for- 
ward matrix elements. Its explicit form is easily 
found by contracting the three-quark matrix cle- 
ment in Eq. (|C.3p with Pq instead of Xa and in- 
serting the resulting expression in Eq. (jC.8p . One 
gets 

^i^2iPa) fdtt{0\e'^''[u'C-f''uq (0)d^(te)|P) = 
4 Jo 



Vx 



Vx 



-^■--^ {^J5N)^ {Px)V2mN 



.^^^JvxJ'dtte-^^■^-^\V^+V,)iJ,N), 



+ ..., (C.12) 

where the ellipses stand for other Lorentz struc- 
-tP-xxs (^^giv)^ Vgrn^v + . .(.C.g) turcs that do not contribute to V^'^^^''^ Inserting 



where [e~*^'^^^(P • x)]. and |^e~'^'^'^-^ ;^] arc the 
leading- twist components for the free fields, defined 
as the solutions of the corresponding homogeneous 
Laplace equation |23- Note that the factor Px in 
the second line in Eq. Ij('.9(l is not included under 
the [. . .]it bracket since (Px) V2 = 1/2(14 -V^z-V^s) 
is a function of momentum fractions only and does 
not contain any dependence on the position vector 
X. To the required 0{x^) accuracy, the leading- 
twist exponents |^e~'^'^^^ ■ ■ ■] it '^^^ ^^ obtained 



everything into Eq. (|C.5ll we arrive at 
(Px) J dx,e-^--P-^V^'^'\x-s) 



(Px) Vxxj / dt 



-iP-xxst 



Vi 



I I Vxx3 / dtte^'*^^^-^(Fi-V^2) 

-^ / Px e-"^^"(-2Fi +V3 + Vi + 2V5) . 



(C.13) 
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In order to solve this equation we expand both sides 
at short distances and obtain the moments of Vj 
with respect to X3 expressed through moments of 
the distribution amphtudes defined as V^ = 

J 'Dxx^^Vi{xi). One finds 



dxsx^V^'^''\xs) = 



1 



n + 1 



v; 



(d)(«+2) 



n + 3 n + 3 

+ ^-(-2^1 + 1/3 + 14 + 2^5)^''^^"+'^ , (C.14) 
n + Z 



up to contributions of multiparton distribution 
amphtudes with extra gluons that have been 
neglected. The corresponding expression for 
Vi i^s) in the momentum fraction space is eas- 
ily obtained by inserting the conformal expansions 
for Vi , . . . , Vg and inverting the moment equation, 
see below. 



The analysis of the u-quark contribution is per- 
formed in a similar way. We consider the matrix 
element 



x" (0|e*^''= [u'{0)C-fa,u^x)] d^(0) \P) = -x°' 



Vi + ^^Vf (")") P, {j,N)^ + V^m^Po. i^j.N)^ 



+ VsTTiAT {jaldN) + V4m%Xa (jsN) + Vem%Xa {M5N) 



and find repeating the same steps that lead to Eq. (|C.5p : 
{0\e^^' [u\0)C^u^{x)] <(0) \P) = (0| [e^^'^ [u\0)C ^iu^ (x)] ^(0)]^^ \P) 



(C.15) 



(C.16) 



— / di 



r^*: r„* 



4 Jo dxadx"^ 



u\{))C ^u' {tx)] 4(0) \P) = (0| [e'J'^ [u\Q)C ^u' {x)] 4(0)]j^ |P) + gluons, 
I 



the only difference being that the term correspond- 
ing to a total translation does not arise in this case. 
For the moments with respect to X2 we get 



da;2X^Vf("'(x2) 



1 



n + 1 



W 



(•u)(n-|-2) 



n + 3 



[Vi - T/2)("^'"+^' 



+ ^-(-2Fi + Vi + Vi + 2y5)(")(«+2) 
n + 2 



(C.17) 

Inserting the conformal expansions for Vi , . . . , Vg 
and inverting the moment equations we find 

1—X2 

Vf(")(x2)= / dXlV^^'{xuX2^-Xl-X2) 



-f(/^C; + AiCr), 



(C.18) 



V^'^''\x3)= I d.Ti^i^(xi,l-xi-a;3,a;3) 
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{jnCJ + X^Cf) 



(C.19) 



with 



CJ=(1 - X2f [113 + 495a:2 - 552a;^ 

- 10A5'(1 - 3x2) 

+ 2^^(113 - 951x2 + 828x2)] , 
Cl=~{l~X2f [13 - 20/f + 3x2 + 10/r(l-3x2)] , 
C^=-(1-X3)[1441 + 505x3 - 3371x^ + 3405x^ 

- 1104x;^ - 24Vf (207 -3x3 -368x3 + 412x^ 

- 138x;^)] - 12(73 - 220^1'') ln[x3]. 



C,^ = -(l-X3) 



11 + 131x3 - 169x2 + 63xi^ 



-30/f(3 + 11x3 - 17x^ + 7x^)] 
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-12(3-10/i^)ln[a:;3] 



(C.20) with 



This result agrees with |35J- 

Similarly we get for the axial-vector functions: 



- n AM{d}, ^ 1 



71+1 



^i' 



(d)(n+2) 



n + 3 

+ ^— (-2A1 - A3 - A4 + 2A5)('')("+2) 
n + 2 



71 + 3 



(yli+A5)('')("+i) 



A/(u) , 



dx2X^Ar^"^(a:2) = 



1 



77+1 



^i' 



(«)(n+2) 



77 + 3 



+ (-2Ai - A3 -Ai + 2A5)(")("+2) 

77 + 2 

which is solved by 



EJ^- 



EI-- 



(1 - x)(1339 + 259a; - 2021a;^ + 1851x^ 
552x'' - 72^5" (1 - xfi?, - 23x) 
24Vf (216 - 99x - 134^^ + 196a;^ - 69a;^)) 
- 12(73 - 220^]^) ln[x] , 

(1 - a;)(53 + 5x - 43a;2 + 21^^ 

- 30(7 ~x-hx^ + 3x^)f^ - 60(1 - x) Vr) 

-12(3-10/{')ln[x], 
^^=31 + 135a: - 138a;2 

- (59 - 483a; + AUx'^){A'l - V{^) , 
i?^=4(l - 3x) - 10(2 - 3a;) (A'' + /{') . (C.25) 

Note that the a;^-corrections do not depend on A2. 
Our results agree with the ones obtained in |45(. 



<^"^(X2) = 



dxiAi (a;i, a;2, 1 — a;i - X2) 



-^(l-:,,f{f^D- + X,Dl) 



24 



Ar''\x3y- 



with 



dxiA^ (a;i, 1 - a^^i - 2^3, 2:3) 
(C.21) 



7:>}^=11 + 45x2 - 2A5'(113 - 951a;2 + 828x^) 

+ 10yi''(l - 30x-2) , 

D5;=29 - 45x2 - 10/i"(7 - 9a;2) - 20/f (5 - 6a;2) . 

(C.22) 

Finally, using the isospin relation in Eq. (2.20) of 
[83j one obtains the T^^ functions in terms of V^^ 
and ^f : 

T/^(«)(a;)=l [F/'W(a;) + V^,*^(")(x) - Af ("Ha;) 

r/'('^)(a;)=l^i^'("'(a;) + Af <"'(x) . (C.23) 

Inserting the above expressions we get: 

T,''^''\x) = ^{f^E- + X,E-), 

^M(.)(^)^ .x^(l-x)^ ^^^^, ^ Aii?,^)(C.24) 



APPENDIX D: ASYMPTOTIC 
DISTRIBUTION AMPLITUDES 

For completeness we present the set of DAs that 
is obtained by setting contributions of higher con- 
formal spin operators to zero, cf. ljB.14p . The sub- 
tlety is that conformal symmetry is broken by nu- 
cleon mass corrections. As a consequence, "kine- 
matic" higher twist corrections of higher spin have 
to be retained in order to satisfy EOM. For the 
relevant parameters we get: for twist-3 



03 = /^: ^3^0, 03+ 



0; 



(D.l) 



for twist-4 

04 = ^ (/iv + Ai) , 01 



V^4 = (/a^ - Ai) , i'i 



1 r 



I ¥^4 






1 
12 



/JV , -04 = 7/. 



N, 



e4°=gA2, ^4^=0, C=0; 
for twist-5 



(D.2) 



6^ = -^(65/^ + 3Ai), ^^ 



1 



[5 In - Ai) , 



1 

18 



^5+ = -T? (65/Ar - 3Ai) , ^|J- = ^ {5/^ + Ai) , 



C5" = ^A2, 4+ = -^A2, 4"--Ja2, (D.3) 
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and for twist-6 
4>l = In, 



-^(5/jv + 6Ai)^ 



(5/jv-2Ai). 



(D.4) 



The corresponding twist-3 DAs are: 

Ti{Xi) = 120xiX2X3fN; 

twist-4: 

V2{X^) = 2X1X2 [5{l + X3)fN + 6X1] , 
A2{xi) = 30xiX2(xi - X2)fN , 
T2{Xi) =4xi2:2A2 , 

V^ixi) = X3 5n + 2xiX2 - 4x3 



(D.5) 



+ 3{xi +X'2+ Xi)j fN - 6(1 - X3)X 

Asixi) = -2(.Ti - a::2)a::3 [5(2 - 3x3)/jv - 3Ai] , 
Taixi) = X3 5(1 + 2a;ia:;2 + 2x3 

-3{xl+ xl+ xl) j /at + (1 - a;3)A2 , 

Tj{xi) ^ X3 5(1 + 2xia;2 + 2x3 

--3{xl+xl+xl)'^fN-il-X3)X2]; (D.6) 
twist- 5; 



V4{xi) = --(28 - 65(1-1 + xl) - 30x1X2 - 13x3 



15x: 



D/n- (l + (a;i-X2)^-x2)Ai, 



Aiix,)^ -(xi-x2)[(37-80x3)/w-6Ai] , 

r4(x,) = -- [2(28 - 65(xi + xl) + 30x1X2 - 43x3 
6 L 



15X3)/jv - (1 + Xi-6X1X2+X2 + 2X3-3X3)A2 



Tsix^) 



2(28 - 65(x? + xl) + 30xiX2 - 43x3 



+ 15X3)/Ar + (1+Xi — 6X1X2+X2+2X3 — 3X3)A2 



V^^Xi) = X3 



^(28-65x3)/iv + 2(l + X3)Ai 



^5(2;,) = ^2:3(1 + X3)A2. 

and twist-6: 

VeixO = ^(5 + 3x3)/Ar - ^(1 - 3x3)Ai 
3 5 

Ae{xi) = --(a;i - X2){5fN - 2Ai) , 



(D.7) 



Teixi) = -(8-6x3)/Ar, 



(D.; 



The corresponding expressions for the x^- 
corrections read: 

2 
V*^^'*' {x) = ^ (l-x) ((215 - 529x + 427x^ 

- 109x^)4- 4 ln[x]Vw + 16(1 -x)^Ai , 
V^^"\x) = ^(1 - x)' [(565 - 417x)/^ - 24Ai] , 
Af"^^\x) = ^(1 - x)^ [(43 + 105x)/^ - 24Ai] , 
r/^W(^) = y (1 - x)3 [(76 - 39x)/jv - 6A1] , 

2 

r/^(")(x) = ^ [((1 - x)(389 - 1051x + 949x2 

- 283x3) + 4 hi[x]) fN + 16(1 - x)3 Ai 1 . (D.9) 



A^ixi) = 2(xi - X2)X3 (5/Ar - Ai) . 



APPENDIX E: QCD SUM RULES 

In this Appendix we update the QCD sum rules 
for the shape parameters of the higher- twist DAs, 
which are defined as |35l 1831 



iuiO)Cj^u{0)) ^^5r{tz'nd){0) \P) = 

^Xiftipz)M^N{P), 

(M"(0)Ca^,M(0)) /75a^''(iz:^d)(0) |P) = 

= X2f.^{pz)M/N{P), 

w(0)C7^75iz D w(0)) ^l^diO) \P) = 

= X,fnpz)M^N{P), (E.l) 



where we have used the notation iz- D~ iz ■ {D — 
D) and for brevity omitted cofor indices. 

The QCD sum rule estimates for f^, /|, /" are 
derived from the consideration of the two-point cor- 
relation functions 



tjd'xe^P^ (0|T{£'^'^ {u^{x)Ca,,,u^{x)) ^,a^^%iz-^d)\x) ^m |0) = ^''^'' ^z^JpZ^^^^^ 
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iJd^xe'P^ (0| Tie'^" {u\x)Cj,u^ (x)) ^^^^{iz^d)\x) i^m |0) = Al^il j|^/J^y + ^'^) + . . . , 
iJd^xe'P- (0| T{£'^'= (^(.t)C7^75*2 D ^(.t))'' 7^rf'(^) ^i(O)} |0) = /"l^il'j^^V ^(^+ M) ^ ^ ^ ^ ^ 



(E.2) 



The dots refer to contributions of excited states and different Lorentz structures that we do not consider. 
We have calculated the correlation functions in ljE.2() using a more consistent factorization approximation 
for the contribution of dimcnsion-8 operators compared to J83l | , which takes into account the contribution 
of nonplanar diagrams. For example, we use 



^.■'c,^'/^'(0|(u'C7m l^c'Sp + ^/3^al ^^')(^''7.g^'')|0)- ^r [195^.^5-/3 - 2(5Maff./3 " 9.<.g^p)] 



216 



kl 



where (uu) is the w-quark condensate and ttiq = (uagGu) /{uu). Following the standard procedure and 
replacing |Aip and IA2P by the corresponding sum rules derived from the diagonal correlation functions 
of the rji and 772 currents, respectively: 



2(27r)4TO^|Aip=exp(TO^/Af2) J M'^Esiso/M^) + -M^Ei{so/M^) + 



4- 



4 TOq 
SAP 



2{2nfm%\^=exp{m%/M') Ul^Esiso/M') + ^M^E,{so/M') 



(E.4) 



we obtain 



^^^M<^E,iso/M^) + lM^E,iso/M^) + ^ (4- f ^ 



/r 






A0E,iso/AP) + ^APEriso/AP) + f 4 - ||^ 



^-AI^E^iso/M^) + lAPE,iso/AP) -' ™' 



10 



2 m^ 
— Q 
3 IP' 



il/SR. 



{so/AP) + |M2£;i(so/M2) + f (4 



3 A/2 



A'^^ 



M/6i?3(so/M2) 



6A/6S3(soA/2) + f M2£;i(so/M2) 



(E.5) 



where 



fc=0 

(E.6) 
In all sum rules AI is the Borcl parameter; we use 
the interval IGeV^ < AP < 2GeV2, with the con- 
tinuum threshold ^/sq ^ 1.5 GeV and also values 
of the condensates at the scale /i^ = 1 GeV^ 

a = -(27r)2(gg)~0.55 GcV^ 



b ^ {2TTf{—G^)c^0A7 GcV 



{qgGq) 



iO.65 GeV^ 



(E.7) 



With these inputs we find the numbers quoted in 
(|B.16ll . These results have smaller errors and su- 
persede the corresponding estimates in Ref. [83 ■ 
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